Curvature-induced quantum behaviour on a helical nanotube 
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We investigate the effect of curvature on the behaviour of a quantum particle bound to move on a 
surface shaped as a helical tube. We derive and discuss the governing Schrodinger equation and the 
corresponding quantum effective potential which is periodic and points to the helical configuration 
as more energetically favorable as compared to the straight tube. The exhibited periodicity also 
leads to energy band structure of pure geometrical origin. 
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Recent developments in nanotechnologyE made it 
possible to grow quasi-two-dimensional surfaces of ar- 
bitrary shape where quantum and curvature effects 
play a major role|2]. Examples include single crys- 
tal NbSe 3 Mobius strips [3], spherical CdSe — ZnS 
core-shell quantum dots [4 , Si nanowire, nanoribbon 
transistors 5. , quantum waveguides [5] and nanotorus[7]. 
Several publications |9j EH EQ E3 EH ED] have 
treated the constrainment of quantum-mechanical par- 
ticles (with applications in, e.g. standard Schodingcr 
equation problems 1 13] and relativistic Dirac equation 
problems Efl El] ) to a two-dimensional surface since the 
original works by Jensen and Koppe, da Costa [TBl[T71ll8j . 
Since two-dimensional systems are an a priori idealiza- 
tion it is reasonable to quantize before constraining the 
particle to the nanotube. As a result a quantum par- 
ticle confined to a two-dimensional surface embedded 
in R 3 experiences a potential that is a function of the 
Mean and the Gauss curvatures of the surface EO H"8] , 
This curvature-induced quantum potential is a geometri- 
cal invariant which property lead the authors [ST] to pose 
the inverse differential geometrical problem: what curved 
surfaces produce prescribed curvature-induced potential. 

Possible physical applications of the above include the 
geometric interaction between defects and curvature in 
thin layers of superfluids, superconductors, and liquid 
crystals deposited on curved surfaces [22 ; the curvature 
of a semiconductor surface determines also an interesting 
mechanism of spin-orbit interaction of electrons [23]; a 
charged quantum particle trapped in a potential of quan- 
tum nature due to bending of an clastically deformablc 
thin tube travels without dissipation like a soliton[24 ; 
the twist of a strip plays a role of a magnetic field and 
is responsible for the appearance of localized states and 
an effective transverse electric field thus reminisce the 
quantum Hall effect [25 . 
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Now let us turn our attention to the geometrical re- 
alization of the helical tube. One can associate with a 
space curve x(s) at any point s along it a moving frame 
consisting of three vectors t-tangent, n-normal and b- 
binormal and evolving along the curve according to the 
Frenet-Serret equations: 

t = a;At, b = cDAb, n = Co A n, (1) 

where Co is the instantaneous angular velocity of the 
Frenet-Serret frame where the arclength s plays the role 
of time. Hereafter the dot denotes derivation with re- 
spect to the natural parameter s. Here re(s) and r(s) are 
the curvature and torsion of the space curve. 

Since Co has a component along t we redefine the frame 
vectors 

N = cos 0(s)n + sin 9(s)b, N = Q A N (2) 
B = — sin 0(s)n + cos 0(s)b, B = flAB. (3) 

We choose 6(s) so that f2 has no component in the direc- 
tion of t. A brief calculation yields 




FIG. 1: The geometry of an infinite helical tube may be 
parametrized by two families of space curves (see equation 
|5| and text). 
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FIG. 2: The cross-section of the nanotube in Fig. [T] 



Now let us mount a disc D rigidly in the reference 
frame where N and B are at rest, i.e. the Fermi- 
Walker frame [27]. The points on the surface may be 
parametrized by 



X(s, <p) = x(s) — po jsin </>B + cos </>N j 



(5) 



The two families of space curves weaving the above 
surface in M 3 are the following. The first is a circle 
parametrized by the angle <fi and is actually the rim of the 
disc that is rigidly mounted to the tangent of x(s) at each 
point in space. The tip of the vector in the disc from the 
central axis to the rim is denoted by po cos tfiN+po sin 0B. 
Its origin coincides with the helical space line x(s). The 
second is given by the lines with tangent passing through 
each point of the first family. Refer to Figure [l] for the 
visual expression of the above construction. 

In this article we will study the properties of the 
Schrodinger equation on that surface. 

The line element is 
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has a dimension of length. 

If we change the parametrization s — > —s and ip - 
this would mean that we evolve the surface backward 
from a certain arbitrary point sq of the infinite space line 
x(s). The torsion r exhibits invariance r(s) = r(— s) and 
the surface element must remain unchanged. 
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h(s) -> h(-s). 



Thus we show that the line element is indeed invariant 

\dX( S ,p)\ 2 = \dX(-s 7 ~p)\ 2 . 

From formulas ^ and ([3 1 we see that at 8(s) — 0, that 
is at s = if so = (see (4|), we have the coincidence 
N = n and B = b. The normal n always points towards 
the axis around which the helix is wound, i.e. it points 





FIG. 3: Density plot of the potentials V e ff(s, ip) on the left 
and V C urv(s,(p) on the right. Along the horizontal axis is the 
natural parameter s, where s € [0, 2n]; <f> is on the vertical 
axis ip G [0,/9o27r]. Here r = K = 1 and po = 10 _1 . Lighter 
regions correspond to higher values of the potential and lower 
probability to find a particle there, respectively. 



inward. From ([7]) it is clear that h(0, 0) = 1 + Pok(0) > 
1 — Pok(0) = h(0,ir). The surface is stretched more on the 
outside thus we have a natural choice of the origin (the 
outer intersection of the ray through n and the cross- 
section of the tube) for the two families of curves (see 
Figure |2). 

Introducing the normal to the surface v from the Gauss 
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we can com- 



pute the linear Weingarten map 
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d s v j ■■ \ <9 S X 
where W is the matrix realizing the map of the tangent 
space in itself 



W 



Po 1 







k(s) cos 0{ s ) + l 



(9) 



With the help of ^ we may compute 



--tr(VF), 



K = kik 2 = det(IF), 



the Mean and the Gauss curvatures of the surface respec- 
tively, where Ki and n 2 are the principal curvatures of the 
surface. They are also the eigenvalues of the Weingarten 
matrix ([£]). Thus we obtain 
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Since we study the resulting Schrodinger equation for 
a particle confined to move on that surface and following 
da Costa an effective potential appears in the Schrodinger 
equation which has the following form: 
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where /i is the effective particle's mass, ft-Plank's con- 
stant; V curv depends on s and p which appear as the 
generalized coordinates on the surface; M = (ki + K<i)j2 
and K = K1K2 are the Mean and the Gauss curvatures 
respectively. For the surface (p>| we obtain 



Here k is the wave number (fc 2 = 2p,E/h 2 ). The cor- 



responding spatial configuration for which equation (16 1 



serves as an effective Schrodinger equation is depicted in 
Figure [T] 

Fixing s — a, straight-farward check provides us with 
the estimate 

V eS (<p = 0) < V eS (tp = ir) 



V cmv (s,(p) = - 



(12) 



From equations ^ and ^ it follows that the surface is 
more stretched on the outside, that is at p — (see Fig- 
ure^, because h(0,0) > h(0,n). The Heisenberg uncer- 
tainty principle states that a particle would have a lower 
energy where the line element is bigger. Our expectation 
is that the probability to find a particle on the outer rim 
of the surface is maximal. This guiding principle will al- 
low us to interpret the appropriate effective Schrodinger 
equation whose potential possesses the above property. 

The Laplace-Beltrami operator (the quantum mechan- 
ical kinetic term) in the coordinate system ^ can be 
written as follows: 
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Here ^> as a solution must be normalized as J \^f\ 2 dS = 1. 
We introduce <& = -\A- V I' and the wave function will be 
normalized with respect to the usual flat norm on a rect- 
angular domain determined by the periodic properties of 
h[e(s),ip], that is f^dtp J**ds\$\ 2 = 1, where s a and 
Sb are such that 9(s a ) = — 6(sb) = 7T. Then: 



and since the probability amplitude follows the behaviour 
of potential the probability to find the particle on the 
outer rim of the surface is greater in accordance with the 
Heienberg's principle. Figure [3] presents the density plots 
of the two potentials. 

Now we will expand the curvature induced effective 
potential V c s and the kinetic operator in series. Next 
we assume k/t ~ 1 and since we want to acquire insight 
into the properties of nanosystems we set pqk = e <C 1 
(po represents the radius of the nanotube, measured in 
nanometers) thus a small parameter naturally arises. 
Now we expand in series the denominators up to first 



order terms in e and equation (16 1 reduces to an effec 



tive two dimensional perturbed Schrodinger equation on 
a rectangular domain in IR 2 : 
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where the perturbing potential is 



fc 2 ff $ = vW$, 
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Now we obtain a differential equation which for the 
helical tube is to be written with k and r constants: 
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Notice, that due to the presence of the squared curvature 
the helical configuration is more energetically favorable as 
compared to the straight tube where this term vanishes. 
This may favor the helical shape in the experimentally 
grown nanotubes |28j. 

It would be interesting to elaborate on the conse- 
quences of the limit e = 0. In this case the geometry goes 
to that of a cylinder and we expect to recover the corre- 
sponding results Indeed the dependence of the wave 
function $ on s is of the form of a standing wave in a one- 
dimensional box stretching to infinity with the boundary 
conditions <!>(s = 0) = <i>(s = L) = as L — ► 00. The 
corresponding eigenencrgies are vanishing 
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What we are left with is a harmonic oscillator equation 
for the dependence of $ on ip. The required periodicity 
of the solution introduces a new quantum number n in 
( 20 1 which lccds to 



E n 



L 



oo. 



This is in exact agreement with the results on the finite 
cylinder [29 . 

Due to the periodicity of the coefficients both in tp and 
s the wave function must be periodic and we may look 
for a solution as a Bloch wave function BUI , that is 
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where S = (2ir) 2 p Q /T is the area of the two-dimensional 
rectangular domain determined by the symmetry of the 
problem and the periodicity of the wave function. In the 
above 



Within the approximation ( 23 ) we obtain 
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Thus the coefficients in the expansion of the Bloch wave 
function (21) have the same order of magnitude as y' 1 ' 



and it suffices to add only few of them in order to obtain 
the behaviour of the wave function on the surface of the 
helical tube. 

It can easily be seen that due to the argument of the 
periodic functions in ( 19 1 we can only have non- vanishing 



contribution along one ray in the inverse lattice, e.g. the 
ray associated with the vector 
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where m s and m v are integers and 
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What is assumed small in this approximation is the per- 
turbing potential of the order of e. It also contains 
derivations which for the reasonably well-behaved wave 



function (21 1 produce no singularities and the order of 



smallness is preserved. 

Inserting this ansatz into (|18[) we obtain 
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k 2 cB \u{K n 
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= V^(K n -K m )u(K m ), 



where are the coefficients of the expansion of the 
perturbing potential (fl9|) along the inverse lattice 
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Thus the direction in the inverse lattice where the sum 



in (21 1 is performed is determined. 



The exhibited singularity at fc, 
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suggests that not only it(0) but also u(K rn ) component 
may be considered as " big enough" . We may write a 
system of equations for the two components directly from 
( 24 1 which gives 



(k+K m y 



u(K m ) ~ V (1 \K m )u(0) = 0, 



V {1 \-K m )u{K m ) + \-k 2 + k 2 cS \ u(0) - 0. (28) 



Equating to zero the determinant of the above sys- 
tem is necessary condition for solvability and produces 



the expressions for the energies due to ( 20 1 in ad- 
jacent zones. Introducing the notation U 2 {K m ) — 
V^{K m )V^\-K m ) we have 



£l,2 = 
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+ U 2 , (29) 



where each root describes an energy band. Here a is 
given by (20 1 and encodes the curvature dependence. It 
is convenient to expand the energy in terms of G which 
measures the difference G = k + \K m in wave vector 
between k and the zone boundary at — \ K m . In the region 



where K 2 G 2 < U 2 
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where £ (±) = a — G 2 ±\U\ , so the energy has 2 roots, one 
lower than the free electron kinetic energy a — G 2 (shifted 
due to the presence of curvature in a) by \U\ and one 
higher by \U\ . Thus the curvature-induced potential has 
created an energy gap 2\U\ at the zone boundary. The 
above expression is valid only when the wave vector is 
very close to the zone boundary. The gap in the en- 
ergy spectrum opened in the transition between the first 
and the second zones in the inverse lattice scales approx- 
imately as 



2\U\ 



(31) 



where k 2 /4 is the curvature-induced potential due to da 
Costa of a space line whose curvature is k. It is of pure 
geometrical origin. In this case the space line is helix and 
the energy gap scales as the ratio between the diameter 
of the nanotube and its radius of curvature, that is e, 
times the energy due to the curved configuration. 



Let us note that the effective mass tensor Mij = 

( W 3 dk"dk ) ^ ^3 IS diagonal. Here the particle ac- 
quires an effective mass due to the interaction with cur- 
vature. This interaction is encoded in the presence of 
the square root containing a in (29 1 which remains af- 



ter twice differentiating with respect to the wave vector 
components in accordance with the conveyed formula for 
the mass. 

In conclusion for a helical nanotube we have obtained 
an effective Schrodinger equation which is periodic. The 
quantum effective potential shows that it is more prob- 
able to find a quantum particle on the outer rim of the 
nanotube. As a whole the helical configuration is more 
energetically favorable as compared to the straight tube. 
The properties of the effective Schrodinger equation are 
discussed within the Bloch ansatz and a gap in the en- 
ergy spectrum is shown to arise in the transition between 
adjacent zones of the inverse lattice. It is geometry de- 
pendent and can be tested also experimentally. 
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